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DEVELOPMENT OF A LASER DOPPLER SYSTEM FOR THE DETECTION 
AND MONITORING OF ATMOSPHERIC DISTURBANCES 

I. INTRODUCTION 


The Marshall Space Flight Center (MSEC) has developed a Scanning 
Laser Doppler Velocimeter System (SLDVS) capable of detecting and monitoring 
the velocity structure of atmospheric disturbances. This system (Fig. 1) has 
been used successfully to detect, monitor, and measure the velocity structure of 
aircraft wake vortices and to provide tracks of the vortex movement. It has also 
been used to obtain vertical wind profiles in a single plane in the atmosphere. 

The SLDVS is a continuous-wave, focused, coaxial, 30.5 cm (12 in.) 

F2 cassegrainian optics system which operates with a CO 2 laser transmitter 
emitting radiation at 10.6 pm. This system is installed in a mobile van and 
consists of a very stable single-frequency CO 2 laser, a modified Mach-Zehnder 
interferometer, a Bragg cell frequency translator, transmit and receive optics, 
an infrared detector, a versatile range and angle scanner, a surface acoustic 
wave frequency analyzer, a signal discriminator, a data algorithm processor, 
displays, recording electronics, and hard copy units. The system is designed 
to have a range coverage from 61 m to 606 m (200 to 2000 ft) and an elevation 
coverage from 3® to 60°. The maximum range and angle scan rates are 7 Hz 
and 1 Hz, respectively. The system detects Doppler velocities up to 61 m/s 
(approximately 200 ft/s) in increments of 0.6 m/s (approximately 1.8 ft/s) 
and provides a line-of-sight velocity spectrum for the range resolution volume 
in space associated with the point at which the system is focused. This velocity 
spectrum is then processed to provide specific velocity/target information as a 
function of position in space and time. 

Aircraft wake vortex and vertical wind profile measurements were per- 
formed with two SLDV systems installed at John F. Kennedy International Air- 
port in New York. These measurements were performed in cooperation with 
the Federal A\iation Agency, the Transportation System Center, and the 
National Aviation Facilities Experiment Center. 



1 1. SYSTEM OPERATIONAL DESCRI PTION 


The SLDVS operational description can best be communicated by tracing 
the coherent output radiation from the CO 2 laser transmitter through the optics, 
to the atmosphere, and back to the data algorithm processor and its displays. 

As shown in Figure 2, the horizontally polarized, collimated 6 mm beam exits 
from the CO 2 laser and impinges on the first beamsplitter which is oriented at 
45° to the beam path. Approximately 10 percent of this beam is passed through 
the beamsplitter to be used as a local oscillator signal. This portion then 
passes through either a Bragg cell frequency translator or a series of attenua- 
tors. In the translated mode, the Bragg cell shifts the frequency by 24 MHz to 
create an offset zero frequency for which positive and negative Doppler shifts 
can be observed. In the nontranslated mode, the beam attenuators reduce the 
local oscillator power to a level for satisfactory opuiaticn of the detector. In 
this mode, positive and negative Doppler shifts and the direction of the velocity 
cannot be determined. The local oscillator beam is then reflected 90° by a 
mirror placed at 45° to the beam patli. The beam then passes through a half- 
wave plate, changing the beam polarization from horizontal to vertical, and a 
recombining beamsplitter where it is mixed with the return signal onto the 
detector. The main portion of the exit beam from the laser is reflected 90° by 
the first beamsplitter and passes through a brewster window that is aligned such 
that it transmits horizontally polarized light and reflects vertically polarized 
light. The exit beam is then reflected 90° and passes through a quarter-wave 
plate that changes the polarization of the beam to right-hand circular. The 
direction of polarization is reversed at the transmit telescope’s secondary 
mirror and again at the primarj' mirror, such that the radiation impinging on 
the scanning flat mirror is again right-hand circularly polarized. Upon reflec- 
tion from the primary, the beam expands to approximately 25.4 cm (10 in.) in 
diameter, strikes the scanning mirror at 45°, and focuses at positions in space 
by movement of the secondary mirror and the scanning flat. The right-hand 
circularly polarized radiation is reflected by atmospheric aerosols in the focal 
volume of the system as left-hand circular and at the same time shifts the 
frequency of the beam by an amount proportional to the velocity of the aerosols 
along the line of sight (Fig. 3). This return signal travels back along the same 
patli as the exit beam until it reaches the quarter-wave plate where, because the 
direction of polarization has been reversed by the atmospheric backscatter, it 
is converted to vertical polarization upon passing through. This vertically 
polarized return signal is then reflected by the brewster window to the recom- 
bining beamsplitter where it is mixed with tlie local oscillator beam. The mixing 
of tliese two beams creates an interference pattern that is imaged by the focus- 
ing lens on the face of the detector. This pattern varies with time according to 
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the difference frequency of the two beams. The detector converts the radia- 
tion energy to electrical energy, and the result is an electrical beat signal that 
is proportional to ihe Doppler shift created by the motion of the atmospheric 
aerosols. This signal is then amplified by the receiver network to the level 
required for operation ox the signal processor. Upon entering the signal pro- 
cessor, the Doppler signal is mixed with a chirp pulse and fed to a surface 
acoustic wave delay line. The serial output of the delay line is a continuous, 
analog spectrum analysis of the Doppler signal where the frequency is propor- 
tional to time. The signal is then sent to a sample-and-hold unit, the output of 
which provides a series of discrete spectrum samples { 104 samples in the 
maximum case). Each sample is representative of a particular Doppler fre- 
quency interval approximately 100 kHz (approximately 1.8 ft/s) wide. Each 
sample is then A/D converted and incoherently integrated by summation aver- 
aging, At the end of a selected integration time (0. 5 ms to G4 ms) , a complete 
Doppler velocity distribution is provided for the aerosols contained in the 
coherent focal volume of the beam for that time and position in space. The 
coherent focal volume or range resolving capability of the focused coaxial optical 
heterodyne systems is expressed in terms of that range element in space that 
contributes a specified portion of the total backscattered signal at various ranges. 
The signal-to-noise ratio of coaxial focused systems that heterodyne the 
scattered return from the atmospheric aerosols is given by 


S/N 


T] /?( 7t) 7T 

2 Bh u 



> 


where 

7 ] = detector quantum efficiency 
= transmitter power 

/3(7t) = backscatter coefficient 
R = radius of the transmitting optics 
B = system bandwidth 
L = particle range 
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h = Planck’s constant 
u = frequency of transmitted radiation 
X = wavelength of transmitted radiation 
f = transmitter/receiver focal length, equivalent to range . 


Integration from Lj to L 2 yields 

- (-¥)] 

The resolution of the system can be defined by calculating the length of the range 
element that produces one-half the coherent backscattered signal received from 
the entire range Lj = 0 to L 2 = ® . Figure 4 gives the one-half power range 
resolution for a 30. 48 cm ( 12 in.) CO 2 laser optical heterodyne system as a 
function of range. This volume is defined as the volume from which one-half 
of the coherent backscattered signal is derived. 

This information, along with the position coordinates, is then recorded 
on high speed, pulsed-code-modulated (PCiM) analog tape at a maximum data 
rate of approximately 2 megabits. The signal next goes through a velocity and 
signal amplitude discrimination process in which diffci’ent velocities and 
amplitudes of the signal may be eliminated, leaving only a selected portion of 
the spectrum (Fig. 5). This information is further examined by the signal 
processor and the following information is extracted: Vmax, the velocity with 
0 highest signal strength; Vpk, the highest velocity with signal above the 
amplitude threshold; Ipk, the magnitude of the highest signal amplitude ( digital 
number between 1 and 256); and N, the total number of velocity cells \vith signal 
amplitude above amplitude threshold. This information is sent to the data 
algorithm processor along with position and time information where it is 
recorded on tape and/or stored in computer memory for processing on a scan- 
frame by scan-frame basis in real time. The results of the data algorithm 
processor are then displayed and automatically hardcopied. 


S/N = 


V P-j, 7 t) X 
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III. VORTEX DATA COLLECTION 


For the John F. Kennedy International Airport aircraft wake vortex test 
program, two SLDVS units were located near the middle marker of runway 31 
ri^t, approximately 121 m (400 ft) on either side of the centerline. This con- 
figuration permitted the two independent sensor units to scan a common area 
perpendicular to the aircraft landing corridor (Fig. 6). This area of prime 
interest for vortex detection and monitoring is indicated in the figure to be 61 m 
(200 ft) either side of the runway centerline and 91 m (300 ft) in altitude. 

In vortex detection, the velocity threshold is set approximately three 
velocity cells, 1. 8 m/s (6 ft/s) , above the peak wind velocity in the scan plane 
and the amplitude threshold is set approximately 6 dB below the system wind 
signal amplitude. This allows only those data points that have high S/N and high 
velocities, i. e. , those that are associated with the vortices, to be sent to the 
data algorithm processor. 

The real time vortex location algorithm processes the raw data to locate 
the vortex centers. The data are screened to ascertain that a sufficient number 
of data points exist to define a vortex center and that the maximum peak velocity 
meets an established minimum. If these criteria are not met, a vortex center 
cannot be defined. When these criteria are satisfied, a correlation region is 
defined that contains all data points within a given radius of the point possessing 
the maximum peak velocity. This region is verified as containing a vortex 
center if B percent of the points possess velocities that are at least A percent 
of the maximum peak velocity where A and B are specified by the operator. If 
this requirement is not met, the point having the maximum peak velocity ( Vpk) 
is rejected as a noise spike, and the entire process is repeated until either a 
noise spike limit is exceeded or a valid region is foimd. When a valid region is 
found, the vortex center is computed by the algorithm in Figure 7. This figure 
also shows the correlation circle drawn about the data point with the maximum 
peak velocity. Once a vortex center is located, the points that were used to 
determine its location are eliminated, and an attempt is made to locate a second 
vortex by repeating the process. After a second vortex center is located, or it 
has been established that one cannot be defined, the chosen output data are 
displayed and recorded. 

The time-based plots ( Fig. 8) display the altitude and horizontal location 
of the port (*) and starboard (® ) vortex centroids as a function of time. When 
only one vortex centroid is found in the scan plane, an ( S) is used to denote its 
position. 



The x-y plots ( Fig. 9) display the vortex centroids in the vertical scan 
plane with scan frame 1 data represented hy an A, frame 2 data by B, etc. , for 
26 frames, at which point the cycle is repeated. The typical time between scan 
frames for the John F. Kennedy International Airjjort operation was 2. 5 s. 

The tabular data outputs (Fig. 10) give tiie fnune number, the number 
of data points contained in that frame, the number of points contained in the 
correlation region for each vortex, the number of noise spikes found while 
determining each vortex center, the minimum and maximum elevation angles at 
which data points were found, the maximum pe:ik velocities in each correlation 
region, the time at which the voitex centers were dcU'cU'd, and die location of 
the vortex centers. 

Scatter plo„s ( Fig, 11) display tlie raw data ))oints in Die jxist-jn'ocessing 
mode, showing tiieir location in die vcrtie.'d scan jilane and Uieir velocities witliin 
3 m/s f 10 ft/s). The vortex centroids arc also displayed. The characters 0 
through 9 correspond to 3 m/s (10 ft/s) velocity increments which start at 0 ni/s 
(20 ft/s) for 0 and go to 36 m/s ( 120 ft/s) for 9. The ten highest velocity data 
points have the actual velocities printed at die Ixittom of the display in feet jier 
second. Those ten velocities arc plotted as characters A Uii ough J. 


IV. WIND PROFILE DATA COLLECTION 

For vertical wind profile measurements, the SLDVS was employed in tlie 
same configuration as in die vortex measurements. Data acquisition was accom- 
plished using the same basic operating parameters. The amplitude and velocity 
thresholds of t'.ic signal processor were set to ensure atmosjiheric wind detection 
and to eliminate signal noise, and the signal integration time was increased 
(from 2 ms to 8 ms) to reduce the total number of data jx)ints in the scan frame. 
The data from the two independent units were Uien jirocessed Ui ensure time 
correlation between scan frames and data points. After a time correlation was 
obtained, the data were examined to find the data points from each sensor that 
could be closely correlated spatially. These points, along with tlicir velocity 
information, were stored and r impared to an input spatird correlation require- 
ment, The spatial correlation parameter requires that data points from the two 
systems must be located spatially at a distance less tJian the spatial correlation 
requirement before tliey can be used to calculate die veitical and horizontal 
velocity components. The j)airs of data jxjinta tJiat met tlie requirement were 
analyzed to provide the vertical and horizontal velocity components associated 
with their common position in the system’s scan plane. Figure 12 shows the 
average horizontal and vertical wind velocity in the '■'•mmon scan area of the 
plane between the two sensors as a function of altitu,;.-. Time correlation, 
spatial correlation, and altitude increments for the c: iculations in this figure 
were 1. 25 s, 2 m (7 ft) , and 6 m ( 20 ft), respectively. 
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V. CONCLUSION 


The Marshall Space Flight Center has developed a Scanning Laser Doppler 
Velocimeter System capable of detecting and monitoring atmospheric disturbances. 
It has been used successfully to detect, monitor, and measure aircraft wake 
vortices of landing aircraft in an airport environment and to measure atmos- 
pheric winds. The data obtained, which include more than 1700 aircraft landing 
operations, have been more than sufficient to permit calculations of vortex 
transport, vortex velocity profiles, vortex circulation strengths and vertical 
atmospheric winds profiles for the scan plane perpendicular to the aircraft 
landing corridor. 
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Figure 1. Scanning Laser Doppler Velocimeter System. 
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Figure 3. Laser Doppler line-of-sight velocity detection. 
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Figure 4. One-half power range resolution versus range for 30. 48 cm (12 in.) CO^ optical heterodyne system. 
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Figure 5. Typical vortex spectrum. 
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Figure 8, Sample time-based plot. 
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Figure 10. Sample tabular mtput. 
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